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Within a mean field plus Random-Phase Approximation formalism, we investigate the collective ex- 
citations of a three component Fermi-Bose mixture of K atoms, magnetically trapped and subjected 
to repulsive s-wave interactions. We analyze both the single-particle excitation and the density 
oscillation spectra created by external multipolar fields, for varying fermion concentrations. The 
formalism and the numerical output are consistent with the Generalized Kohn Theorem for the 
whole multispecies system. The calculations give rise to fragmented density excitation spectra of 
the fermion sample and illustrate the role of the mutual interaction in the observed deviations of 
the bosonic spectra with respect to Stringari's rule. 
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I. INTRODUCTION 

Shortly after the achievement of Bose-Einstein con- 
densation of dilute, magnetically confined alkali atoms, 
an experimental quest began, seeking to trap and cool 
fcrmionic isotopes of lithium and potassium below their 
Fermi degeneracy temperatures Tp The first re- 

alization of a weakly interacting, degenerate Fermi gas 
of K atoms was reported by DeMarco and Jin ||, and 
most recently, a mixture of bosonic and fermionic Li iso- 
topes has been sympathetically cooled to about 80% of 
the corresponding Tp ||. On the theoretical side, vari- 
ous studies of the structure and equilibrium density pro- 
files of harmonically confined fcrmions were anticipated 
in the last years JtJ , as well as thermodynamic properties 
||. The coexistence of trapped bosonic and fermionic 
isotopes of a given element received particular attention 
|5] |l5| , mostly seeking to establish the effects of quantum 
degeneracy and mutual interaction on the respective den- 
sity profiles, energetics and stability properties. 

Since the early days of Bose-Einstein condensates, the- 
oretical and experimental investigation of the excitation 
spectra of these objects constituted a field of active re- 
search (see, for example, Refs. |l6|,[l7j for substantial re- 
views) which largely enriched with reciprocal feedback. 
In spite of the lack of experimental data concerning 
trapped degenerate fermion gases, theoretical study of 
their collective oscillations also developed according to 
two viewpoints. On the one hand, a hydrodynamic ap- 
proach derived from the balance equations for the mo- 
ments of the Wigner distribution function permitted to 
compute the sound mode spectrum of one Q and two 
p9[ hyperfine fermion species. Such an approach, which 
can be generalized for anisotropically trapped systems 
pp| and was applied as well to predict the characteristics 
of scissor modes in a superfluid Fermi gas , conceals 
the hypothesis of underlying local equilibrium of the os- 
cillating gas, i.e., the validity of the Fermi gas equation 
of state. On the other hand, since collision rates in a very 
cold and strongly dilute gas are small for moderate num- 



bers of confined particles |P|^2|] , one may expect that a 
zero sound regime shows up at low excitation energies 
and temperatures in these systems; in fact, collective 
energies of a single-species fermion gas were estimated 
within a sum-rule approach p3| for the lowest multipo- 
larities. Recently, we have proposed a Hartree-Fock (HF) 
plus Random-Phase-Approximation (RPA) description 
of multipolar excitations of two hyperfine components 
of 40 K atoms, mutually coupled through a repulsive s- 
wave interaction (24) . The procedure proved adequate to 
obtain the trend of collective energies and transition den- 
sities for multipolarities L = 0, 1, 2, both for the density 
and spin density oscillations, and the results were satis- 
factorily checked against a simple model, which permits 
analytical predictions in the case of identical trapping 
potentials and particle numbers for either component. 
Moreover, the philosophy turned out to be equivalent to 
the one undertaken by Bruun pSfl for Li atoms interact- 
ing through an attractive coupling. 

Due to their coexistence in natural samples, the possi- 
bility of simultaneous trapping and thermal equilibration 
of both bosonic and fermionic isotopes of a given alkali 
within their respective degeneracy regimes, makes room 
to an interesting field for theoretical anticipation, where 
the study of excitations of these binary mixtures consti- 
tute an open issue. This subject has not been pursued in 
depth, as compared with the interest deposited in struc- 
tural characteristics and energetics. A three fluid model 
for trapped fermions coexisting with both condensed and 
noncondensed bosons was applied in Ref. jl(| to the study 
of the temperature dependent profiles, and in Ref. 
to outline a study of density fluctuations within the hy- 
drodynamic regime. Some features of collisionless, zero 
sound oscillations in the Landau limit have been dis- 
cussed in Ref. ]l4| and in Ref. [^6) , a sum-rule approach 
was applied to extract the lowest collective zero sound 
energies of an isotopic mixture of K atoms. 

The purpose of the present work is to contribute to 
this wide field reporting calculations of the zero sound 
spectrum of trapped, weakly interacting K isotopes - ei- 
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ther bosonic 39 K or 41 K, together with two hyperfine 
species of 40 K. The theoretical frame is the mean field- 
plus-RPA for the coupled threefold system at zero tem- 
perature, which is exposed in Sec. [o] within a simplified 
scope. The results here presented then indicate the trend 
to be expected at very low temperatures, when quantum 
degeneracy overcomes the effects of interactions in the 
fermion sample, and should be revised if the fermionic 
gas lies in its superfluid regime. The procedure can be 
extended to finite temperatures, at some computational 
cost. Specific calculations are analyzed in Sec. Ill, while 



the conclusions and overview of this paper are the subject 
of Sec. EvL 



II. THE PHYSICAL SYSTEM 

As anticipated in the Introduction, in this work we 
study a cold dilute K gas with three components, con- 
fined in a spherical harmonic potential. Two of the com- 
ponents, hereafter called species 1 and 2, are fermions 
with different spin projections and the third one is a bo- 
son B species. Although in principle, different compo- 
nents subjected to the same magnetic field would effec- 
tively feel different trapping potentials, we have assumed, 
for simplicity, that all atoms experience the same confine- 
ment 
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where to is the atom mass, in other words, we neglect mi- 
nor mass differences among potassium isotopes, as well as 
variations in the trap frequency due to the distinct spin 
projections. In addition to the external field, the atoms 
can interact, but due to the diluteness and temperature 
of the gas sample it is enough to consider only s-wave 
interactions ]17| . Since our gas includes two types of 
identical fermions, not all s-scattering processes are per- 
mitted; in fact, each fermion species is a free Fermi gas 
which can only interact with atoms of a different kind. 

Taking into account that the s-wave interaction can 
be represented as a contact potential V s — gS(r — r') 
with g = 4ira s H 2 /to, being a s the scattering length of 
the process, the mean-field Hamiltonians for the three 
species can be written as 
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being gs, gi, g 2l and <?/ the interaction strength between 
bosons, species 1 and B, species 2 and B, and species 1 
and 2 respectively. 

These Hamiltonians provide the single particle wave- 
functions and energies necessary to find the equilibrium 
densities pB,Pi,p 2 , with ps = Nb |</>o| 2 and 
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for a = 1, 2, being N(e) the Fermi occupation numbers 
at given temperature T, and Efo the Fermi energy of 
species ct, satisfying 



(2.4) 



it / 



The particle eigenspectrum is the solution of the cou- 
pled equations 



H^r, 



l<Pnlm, 
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(2.5a) 
(2.5b) 



It is worth mentioning that Eq. (2.5a) for n = 0, 1 = 
is the Gross-Pitaevskii equation (GPE) for Nb atoms 
with chemical potential p — eo in an external potential 
V = Vtrap + + .92P2)- The higher modes are the 

excitation modes of the GPE which must not be confused 
with the Hartree-Fock-Bogoliubov (HFB) modes. 



A. Collisionless approximation 

To obtain the spectrum of collective excitations of the 
compound system we resort to a random phase approxi- 
mation, which treats the three species on the same foot- 
ing. In a previous work |24|] , we have derived two equiv- 
alent sets of two equations each, for the density fluctu- 
ations 8p aa , being a, a' either 1 or 2. These represent 
the response of the density p a to an initial perturbation 
of p a > . The diagonal and off-diagonal density fluctua- 
tions can be later mapped onto spin symmetric and spin 
antisymmetric ones for each species and for the whole 
gas. A similar approach in the presence of the bosonic 
component would lead to nine equations which include 
transition densities Sp BB ,Sp aB ,8p B<T and corresponding 
ones for the other fermion species. In this work, we shall 
consider identical fermion numbers Nfo in each species, 
and we are mostly interested in total transitions densi- 
ties, namely 



s Pa = s P aa + s P aa ' +s P aB 1 

Spb = Sp BB + 8p Ba + 5 P Ba ' 



(2.6) 



Accordingly, from the full system of equations we can de- 
rive a simplified set, which turns out to be similar to the 
one developed by Minguzzi and Tosi for Bose-Einstein 
condensates at finite temperatures [0. The details are 
outlined in Appendix A. In this frame, a density fluc- 
tuation evolves in a perturbed H F po tential according 
to the linearized Hamiltonians va (2/2). In terms of the 
particle-hole (ph) propagators for each spin species in the 
HF approximation, Gg , we get for the transition densities 
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6p a (r) = Sp° a (r) 

+ J GZ(r,r')[ 9<T 5p B (r')+g f 6p a ,(r')]d 3 r', 
Sp B (v)=Sp B (r) 



+ / G s (r,r') 



g B 5p B (r') +^g a Spa(r') 



(2.7) 

where Sp® is the density fluctuation induced by the ex- 
ternal probe in the stationary HF mean-field approx- 
imation 



5p°(r) = / 0\v')Gl(vy)d 3 r' 



(2.8) 



and the propagators read 

^(rO^r)^)^!-') 
il — (ej — p)/h + if] 



G* = N B J2 



fo(r)05(r')^-(r')^(r) 
+ (ej — p)/fr + it] 

N{e a k ~ e F a) - N(e? - e Fo 



(2.9a) 



(2.9b) 



where each single particle label stands for n,l,m, and 
the condensate mode is excluded from the summation in 



fl2.9aD . 

It is clear that in the present situation, i.e., Nfi — Np2 
and identical trapping field for both hyperfine compo- 
nents, one has a single fermion spectrum ef u = ef tl , 
one ph propagator Gq = Gq and same fluctuations 
8 p a = 5pF- Due to the spherical symmetry imposed 
to the tra ppin g potential and for multipolar operators 
O', Eqs. (2.7) can be rewritten in terms of multipolar 
components as 



8p l a (r) = Spl l (r) + 
[da Sp l B (r 
Sp l B (r) = Sp B l (r) + 
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where the expressions for the Z-polar propagators G l ; are 
given in Appendix B and the total density fluctuation for 
type i-atoms reads 



5 Pi (r)=^y i0 (f)^(r 



(2.11) 



The whole process consists of three stages. First, one 
has to find the equilibrium densities, which means to 
solve the GPE coupled to the fermionic system of HF 
equations. This calculation also yields the quasi-particlc 
spectrum of the fermions. Then, one computes the GPE 
excitations for t he bosons, and finally, one solves the RPA 
equations ( 2. 10] ) for the density fluctuations with given 
I. After computing the transition densities, we calculate 
the total susceptibility (x) and dynamic structure factor 
(S) of the system as 



x l m = XB +x[+xi 

S l (n) = --lm [x l ], 



(2.12) 



with x\ the susceptibility of each component in the in- 
teracting system, 



r 2+l 5p\{r)dr 1^0 
r 4 <5p°(r)dr I = 



(2.13) 



From the poles of the real part of x l or corresponding 
peaks in S we can extract the excitation energies of the 
collective modes. 



III. RESULTS 

In this work, we analyze the spectrum of elementary 
and collective excitations for varying concentrations of 
the fermionic species. For this sake, we tune the param- 
eters close to the experimental setup of DeMarco et al. 

i.e , we choose to — 2n x 70 s _1 , a\ = a,2 = 2.519 
nm, af = 8.308 nm and a B = 4.292 nm, corresponding 
to 40 K (fermions) atoms and 39 K (bosons). 

Th e co mputation of the single-particle (sp) spectrum, 
Eq. (2.5), was carried out by expanding their eigenso- 
lutions in a basis of harmonic oscillator functions with 
1800 radial wavefunctions. In fact, due to the spherical 
symmetry of the trap, we have just diagonalized several 
small systems for different I (I < 60). Convergence was 
achieved with an accuracy of at least 10~ 10 in the ener- 
gies, and the coincidence between p and eo was checked 
to about 10 -5 . Then, the numerical solution of Eq. ( 2.10| ) 
for each and L was performed by discretizing the prob- 
lem in a spatial mesh containing up to 250 points and 
extending up to R « 30 pm (approximately fifteen oscil- 
lator sizes); the resulting linear system was solved with 
standard routines from the LAPACK Fortran library. 

In the next subsections we analyze both the quasi- 
particle spectrum and the selfconsistent RPA collective 
modes. We use the notation Np for either Np a and g a 
for the common coupling g\ = gi- 



being Yi m standard spherical harmonic functions. 



3 



A. Hartree-Fock and GPE excitations 

Within the mean field approach to the coupled system 
we cal cula ted the quasi-particle excitations according to 
Eqs. (2.5). The main features of the HF fermionic en- 
ergies are the following: for low numbers of fermions in 
the trap, Np < 10 4 , the excitations energies of the pure 
fcrmion system are basically the harmonic oscillator ones. 
When we turn on the coupling to the bosonic component, 
the repulsion with the condensate increases the Fermi en- 
ergy and lowers the slope of the dispersion relation. The 
change in slope is more significant, the larger the dif- 
ference Nb - Np; indeed, it becomes unimportant when 
these numbers are comparable or slightly negative. 

As shown by M0lmer |J, for high number of bosons, 
their kinetic and interaction energies with the fermions 
can be neglected in Hb, yielding a Thomas- Fermi ap- 
proximated (TFA) |28| bosonic particle density 

PB = N B \M 2 = —(M ~ VWp)e(// - 14rap), (3.1) 
9B 

where Q(x) is the step function. This is translated, aft er 
replacement in the fermionic Hamiltonian [cf. Eqs. (2.2)], 
into a wider harmonic trap, shifted upwards according to 

V c ff = Vti a p + Vtrap)©(A* ~ Vtrap), 

9b 

— M + ) ^r<R B 

9B 9B 
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otherwise 



(3.2) 



with the Bose radius in the TFA, R B = \/2/i/mw 2 . 
For the values of the coupling constants here consid- 
ered, the strength of the restoring force experienced by 
fermions inside the bosonic cloud is reduced by a fac- 
tor (1 — ga/gs) ~ 0.41. Consequently, the Fermi radius 
Rf = \J1sFjnvJ 1 acquires an effective value 



R F 
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2(e F - ng a /g B ) 
(1 - 9*/9B)muj 2 



(3.3) 



which increases with respect to its size in a pure 
fermion system. The low-lying quasi-particle excitations 
— i.e. those with energies below the chemical potential 
(X of the bosons — are then expected to be sensitive to 
the renormalized trapping strength, whereas the higher 
energy ones, which extend further beyond the central re- 
gion, are less affected by the effective trapping. In Fig. [j] 
we plot the dispersion relation for Np — 9880, and sev- 
eral values of N B - The corresponding values of the bo- 
son chemical potential are indicated in the caption. We 
observe that, in agreement with the effective trapping 
picture, when the bosonic cloud is large enough, the low 
lying modes increase with a smaller slope than the modes 
with energies higher that fi. Although the fermionic spec- 
trum is altered by the interaction, the bosonic counter- 
part is almost unchanged; indeed, in addition to the dif- 
ferent number of particles in each isotope, the fermionic 



character of the 40 K atoms gives rise to quasi-particle 
excitations rather sensitive to variations in the effective 
trap. The bosonic spectrum is very different; it actually 
consists of the GPE single-particle excitations which will 
build the collective energies in a RPA calculation. 



B. Random-Phase Approximation 



The general aspect of the response (2.12) can be in- 
terpreted from its constituents. Without the boson- 
fermion coupling, i.e. g\ = #2 = 0, the bosonic re- 
sponse consists of the well-known collective excitations 
enclo sed i n the zero temperature Gross-Pitaevskii equa- 
tion (2.5a), namel y the HFB spectrum mentioned imme- 
diately above Sec. [I A . For high number of bosons, this 
is well reproduced by Stringari's hydrodynamic formula- 
tion ps| , which provides an analytic expression for the 
multipolar excitation energies. Although analytic for- 
mulae have been also found for systems of interacting 
fermions p5| , pl[ |, they hold in the hydrodynamic regime 
and for population ranges far from those we address here, 
so we compare our results with earlier numerical calcu- 
lations carried within the collisionless regime. 

We have analyzed the dynamic structure factors and 
fluctuations created by multipolar external fields with 
L = 0, 1, 2. In Fig. ^ we show the monopolar spectrum 
for several values of Nb and Np. The essential features 
of the monopolar response can be summarized as follows. 

1. Fixed Np : for high Nb (higher rows), the spec- 
trum is dominated in amplitude by the bosonic col- 
lective modes with frequencies fi « \/5, Vl4,...[cf. 
Ref. |§8)]. However, strong fragmentation occurs 
near the fermionic RPA poles, which can be visu- 
alized in the lowest row. As we lower Nb, this 
fragmentation is reduced and two sets of two peaks 
each is distinguished, each set corresponding to the 
bosonic and fermionic constituents. The bosonic 
peak separates from the Stringari spectrum towards 
the ideal gas result Q = 2nui, revealing the unim- 
portant role of both the boson-boson and boson- 
fermion interaction in this very dilute regime. The 
fragmentation in the fermionic part cannot be seen 
due to its weakness and the current energy scale. 

2. Fixed Nb- one realizes that increasing Np not al- 
ways yields a more fragmented spectrum. In fact, 
some peaks appear better defined for larger val- 
ues of Np, while other states seem to spread over a 
larger energy range. On the other hand, one should 
keep in mind that the scale in this plot is governed 
by the boson strength. If this scale is released, 
one may observe the pure fermion RPA peaks (i.e., 
those with Nb = 0) within restricted intervals, with 
fragmented structures for sufficiently high Np. 

The results for the dipolar mode [L = 1) displayed 
in Fig. U are somewhat different. As we can see, for all 
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values of Nb and Np the highest peak lies at = ui, 
corresponding to a center-of-mass (Kohn mode) oscilla- 
tion of the whole system. It is known |2^,|2^,^| that both 
the pure boson and pure fermion systems possess a well 
defined Kohn mode within the RPA formulation, which 
is recovered in the fully coupled calculation, as can be 
seen in Fig. ^. In fact, in the numerical calculations the 
precise location of this pole has been used to check the 
computational code. It should be mentioned, however, 
that in a actual experimental setup, the slight differences 
between the trap forces for each constituents would dis- 
tort these modes, giving rise to coupled oscillation of the 
components. In addition, the dipolar spectrum is highly 
fragmented around fi = 3a;, which coincides with an os- 
cillator mode for 1 = 1. As in the monopole case, the 
important fragmentation seems to disappear as we lower 
either Nb or Np. 

The quadrupolar mode exhibits the strongest fragmen- 
tation. Even for low values of Nb and Np the spectrum 
presents several major peaks whose number grows at 
higher densities. At the highest population of all species, 
the largest peak corresponds to a boson excitation with 
in agreement with the sum-rule estimate of 
Ref. p9] . As we diminish the number of bosons the peak 
is shifted upwards towards the ideal gas value f2 = 2ui. 
Another feature not seen in other multipolarities is the 
existence of peaks with substantial strength at rather low 
energies near w 0.3a;. 

The shape of the transition densities deserves its own 
study. In Fig. ^| we depict the monopolar fluctuations cor- 
responding to the frequencies where the maximun transi- 
tion probability is attained, each curve has been indepen- 
dently scaled to fit in the plot. Apart from differences in 
scale, the shapes of the bosonic fluctuations are very sim- 
ilar for all Nb and Np, displaying in any case a volume- 
like excitation with nonvanishing density at the origin. 
However, two kinds of fermionic excitations can be seen, 
a volume-like one, similar to the bosonic mode, and an- 
other which oscillates between the boundaries i?,Feff and 
Rb of the isotopic particle densities. These two shapes 
can be correlated to the spectrum [see Fig. ^ : while the 
volume-like fluctuations occur for essentially pure boson 
peaks and in configurations where the bosonic cloud sur- 
rounds the fermionic component, the surface-like mode 
appears near the pure fermion excitation. The fermionic 
atoms extend outwards with respect to the bosonic den- 
sity and the excitation bounces between the fermionic 
and bosonic density boundaries. To support this idea we 
have tabulated the radius of both densities in Table QL 

The monopolar response contains other types of den- 
sity fluctuations, some of them correspond to important 
poles in the response; one can observe not only volume- 
like bosonic excitations, but also surface-like ones, in spa- 
tial counterphase to the fermionic ones for the frequencies 
shown. 

As expected from the Generalized Kohn Theorem, the 
dipolar fluctuations at trap frequency do exist and cor- 
respond to a small displacement of the density yielding 



a density fluctuation proportional to the gradient of the 
density, i.e. 

5pr i (r,n = a;)oc^. (3.4) 

For high enough Nb, according to the TFA, the 
bosonic density shape inside the condensate is parabolic, 
giving a linear bp. This is seen in Fig. ^| where we plot 
the fluctuation densities corresponding to these modes 
for Ng = 10 6 and Np = 969. The shape of the fermionic 
fluctuation is not as simple, and in fact, should not have 
a linear dependency with position since the TFA den- 
sity profile carries a power 2/3. Moreover, the oscillatory 
shape shown in the plot agrees very well with the deriva- 
tive of the fermionic density calculated within the HF 
approximation. 

Finally, the behavior of the L — 2 density fluctuations 
is similar to the monopolar ones in the sense that for large 
number of bosons and few fermions, the highest peak in 
the response corresponds to a fluctuation which extends 
over the bulk of the cloud; whereas for low boson num- 
ber and/or many fermions, the fermionic fluctuation is 
mainly surface-like. The difference with L = lies in the 
fact that the quadrupolar oscillation vanishes smoothly 
at origin. 

IV. CONCLUSIONS AND OVERVIEW 

In this work we have designed a mean field plus RPA 
formalism for a trapped Fcrmi-Bose mixture, i.e., a three 
component gas at vanishing temperature. We have ana- 
lyzed both the single-particle excitation spectra and the 
collisionless collective modes. On the one hand, the 
fermionic single-particle spectra can be interpreted in 
terms of an effective confining field induced by the pres- 
ence of the bosonic atoms plus the trapping potential 
|p[-|rT|. This effective field produces a packing of sp in 
the lowest branch of the spectrum. On the other hand, 
the bosonic part of the spectrum is mainly unaffected by 
a moderate fermionic cloud. 

The poles of the monopolar and quadrupolar response 
are shifted away from either the pure fermion or pure 
boson results and the whole spectrum is highly frag- 
mented, specially for the quadrupolar probe. In addi- 
tion, for the dipolar excitation, the formalism proves to 
be theoretically and numerically consistent with the Gen- 
eralized Kohn Theorem and the density fluctuations are 
proportional to the derivative of the corresponding equi- 
librium densities. The shape of the monopolar fluctua- 
tions reflects two kind of behaviors. For high number of 
bosons, both the fermionic and bosonic fluctuations are 
volume-like, extending up to their respective densities 
edges, while for similar number of particles, although the 
bosonic fluctuation remains volume-like, the fermionic 
counterpart bounces between both density borders. 
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Finally, this formalism can be extended to finite tem- 
peratures where the system may be found partially de- 
generate, i.e. with the bosonic and/or fermionic compo- 
nents out of the quantum degeneracy regime. This new 
regime could give rise to an interesting and rich field of 
research. 
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Under the conditions specified in this work, in other 
words, identical trapping potentials and same number of 
fermion of each hyperfine type, the fermion propagators 
are identical and one can easily construct the equations 
for the total density fluctuations of the fermionic and 
bosonic atoms, summing over the hole partner 
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(A5) 



The final expressions are displayed in Eqs. (2.7). 



APPENDIX A: RANDOM PHASE 
APPROXIMATION 

We assume that the trapped, three component atom 
system is excited by the action of an external field, so 
that ph pairs (a(3) involving, in principle, both hyperfine 
fermion species and the bosons, are created. We also sup- 
pose that a ph effective interaction V^^ S , that permits 
propagation of elementary excitations carrying energy f2 
to all orders in the interaction, gives rise to the dressed 
RPA propagator for ph pairs |31[ 



(Al) 
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We consider longitudinal excitations involving prop- 
agation of ph pairs of the same spin kind, created by 
multipolar operators of the form 



r L Y L0 (9,cp) L^O 



r 4 r 00 



L = 0, 



(A2) 



and an effective interaction which couples ph pairs of i) 
different fermion species, expressed as = V£™ a \ 

ii) each fermion species and the bosons, V° B = V^ BB 
and iii) bosons, V^ B . The free ph propagators involved 
in longitudinal densi ty fluctuations are then diagonal in 
spin space and Eq. (Al) gives rise to nine coupled equa- 
tions. After convoluting the propagators with the exter- 
nal field, one gets an equivalent system for the density 
fluctuations, which in matrix form can be cast as namely 

Sp" = 6p%° + G%° (g f 5 P °'° + g a 5p B °), 

5p°'° = G£ V (g f 5p™+g a 6p B °), 

b P B ° = G BB ( 9<7 5p™ + g a , 5p°'° + g B 6p B °), 

5 P ° B = (g f h°' B + 9o 5 P BB ), (A3) 

plus four similar equations corresponding to interchang- 
ing a and cr', in addition to the modified boson RPA 
equation 



APPENDIX B: MULTIPOLAR DECOMPOSITION 

To evaluate the free ph propagator we us e th e 3D ba- 
sis functions calculated according to Eq. (2.5). These 
eigenfunctions can be written, for each species, as 



(f>nlm(r) = Rnl(r) Yl m (f), 



(Bl) 



where Yi m are spherical harmonic functions and R} n l the 
corresponding radial eigenfunctions. 



Decomposing the propagators Eqs. (2.9) as 
Q(r,r',^) = ^Q i P i (f.;'), 



(B2) 



with Pl{u) Legendre polynomials of order L. The mul- 
tipolar component G l Q L read 

GZl = J^y E R n I W R n' V I (O K' V </) 



being 
and 



(2/ + l)(2l> + l)\(l0l>0\L0)\ 2 x a ni, n 'i>m, (B3) 
N(e n i - e F a) - N(s n > v - e F<J ) 



(B4) 



u - (enL - p)/h + ir) u; + (e n L - p)/h + ill ) ' 

(B5) 



where the primed summation in Eq. (B5) indicate to ex- 
clude the condensate mode, i.e n > 1 for L = 0. Simi- 
larly, the /-polar density fluctuation read 



4tt 



2L+ 1 




12+L rii 



Gl L {r')dr' L^O 
r' 4 G l 00 (r')dr' L = 0. 



(B6) 



G 



APPENDIX C: SIMPLE PROOF OF THE KOHN 
MODE 

The Kohn theorem states that, independently of the in- 
teractions, a many-body system subjected to a harmonic 
confinement possesses an excited state corresponding to 
an oscillation of its center-of-mass at the bare oscilla- 
tor frequency = cj . Consequently, the fulfillment of 
Kohn's theorem is a desired feature of any approximate 
formulation. 



The RPA equations (A3) and (A4) can be cast in a 
9x9 matrix form, namely 



e -Sp = Sp Q , 
with Sp = {5 P n , Sp 21 , 5 P B \ Sp 12 , Sp 22 , 5p B2 , 




(CI) 
■)*, and 



(C2) 



being 



/ 1 

~9fG% 

V -9i G B 



-9fGh 



-92 Gi 



-9i Gl \ 
-92 G 2 
1 - 9b G B ) 



(C3) 



The collective states of the system are located at the 
frequencies f2 which give a singular matrix e, implying a 
singular A. In addition, finding an eigenmode of A with 
zero eigenvalue we prove that A is singular. More over , 
solving this eigensystem is equivalent to solve Eq. (2.7) 
with 8p® set to zero. 

We propose the solution 



r._, X / Upi Up 2 



dpi dp 2 dp B 



dxj ' dxj ' dxj 



(C4) 



(for any direction j) and using the fact that G l Q is built 
of single-particle excitations corresponding to the mean- 
field Hamiltonians ([2.2]) we find [B2] 



/ 



Jg5(t,t',u> 
G B {r,r',u) 



9a 



df>E 
dx>. 



9f 



9pa 



3„/ 



9b 



dpB 
dx' 4 



Edpa 

a J 



d 6 r' = 



Opa_ 

dxj 

dps 

dx-j 



(C5) 



which correspond to the requested eigensolution of 
5p = 0. 



Also at the Consejo Nacional de Investigaciones 
Cienti'ficas y Tecnicas, Argentina. 
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FIG. 2. Monopolar dynamic structure factor in the RPA 
approximation (in log 

scale and arbitrary units). Each row corresponds to, from 
top to bottom, N B = 10 6 , 10 5 , 10 4 , 10 3 ,0; while each column, 
from left to right, N F = 969, 2925, 9880. The vertical dotted 
lines show Stringari's energies psj for the bosonic system. 




FIG. 1. HF fermionic excitation energies (in units of hu) 
for Np = 9880 and I = 0. The different symbols indicate dif- 
ferent number of bosonic atoms: (o): TVs = 10 6 , fi = 32.5hw; 
(>):N B = 10 5 ,/x = I3.2hu; (<): N B = 10 4 ,^ = 5.7hto; and 
(*):Nb = 10 3 , fj, = 2.9hu. The dashed-dotted line shows the 
result for g 1 = g 2 = 0. 
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FIG. 3. Idem Fig. H for the dipolar response 
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FIG. 4. Idem Fig. ^ for the L — 2 structure factor. 
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FIG. 5. Transition densities for the L — response (in ar- 
bitrary units) . The plots show the real parts of the transition 
densities at the highest peaks in Fig. |^ for each (Nf, Nb) 
pair. The solid lines show the fermionic fluctuation while the 
dashed ones indicate the bosonic part. Each plot corresponds 
to the same ordering as in Fig. ^| The positions are expressed 
in terms of the oscillator size d = ^Jh/rnuj » 1.9/im. 




5 10 

r/d 



FIG. 6. Dipolar fluctuation at the oscillator frequency for 
TVs = 10 6 and Nf — 969 (in arbitrary units). Solid and 
dashed lines correspond to the fermionic and bosonic fluctua- 
tions respectively. Each curve has been scaled independently. 

TABLE I. Parameters of the stationary densities. Rb, 
RfcH, ^ and ef respectively correspond to the radius of the 
bosonic and fermionic densities and the chemical potential of 
the bosons and Fermi energy. 



Nb 


Rb /d 


Rf off /d 


fi/hui 


£F /hbO 


N f = 969 










10 6 


8.3 


7.0 


32.5 


30.5 


10 5 


5.2 


6.0 


13.0 


20.3 


10 4 


3.8 


6.0 


5.4 


18.5 


10 3 


3.0 


5.8 


2.6 


18.5 


N F = 2925 










10 6 


8.4 


8.4 


32.4 


36.1 


10 5 


5.3 


7.2 


13.1 


27.6 


10 4 


4.0 


7.0 


5.5 


26.4 


10 3 


2.5 


7.0 


2.7 


26.4 


N F = 9880 










10 6 


8.0 


9.2 


32.5 


46.0 


10 5 


5.3 


8.7 


13.2 


40.1 


10 4 


3.8 


8.5 


5.7 


39.3 


10 3 


2.4 


8.5 


2.9 


39.3 







